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ABSTRACT 

We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for 
galaxies with stellar masses as low as 10 5 M 0 by combining star formation histories of 
Milky Way satellite galaxies derived from deep Hubble Space Telescope observations 
with merger trees from the ELVIS suite of N-body simulations. This approach extends 
our understanding more than two orders of magnitude lower in stellar mass than is 
currently possible by direct imaging. We find the faint end slopes of the mass functions 
to be a = 1.42^q q 5 at z = 2 and a = 1.57^06 at z = 5, and show the slope only 
weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of 
direct detection studies that suggest slopes as steep as a = -1.9 at these epochs. Such 
a steep slope would result in an order of magnitude too many luminous Milky Way 
satellites in a mass regime that is observationally complete (M* > 2x 1O 5 M 0 at 2 = 0). 
The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end 
slopes that are consistent with our results, but with a lower degree of precision. This 
work illustrates the strong connections between low and high-z observations when 
viewed through the lens of ACDM numerical simulations. 

Key words: cosmology: theory - galaxies: dwarf - galaxies: high-redshift - galaxies: 
evolution - galaxies: luminosity function, mass function, galaxies: Local Group 


1 INTRODUCTION 

The evolution of the galaxy stellar mass function, from high- 
redshift to low, provides an important benchmark for self- 
consistent galaxy formation models set within the Lambda 
Cold Dark Matter (ACDM) framework (White fe Rees|1978 


Somerville fe Dave||2014| |Vogelsberger et al.||2014| |Schaye] 


et al. 


20151. In particular, the low-mass slope of the stel¬ 


lar mass function (hereafter, afaint) is a crucial gauge of 
how and when feedback becomes important in suppressing 
galaxy formation in small dark matter halos. Indeed, many 
of the best models today predict that the galaxy stellar mass 
function should be much steeper at high redshift than low 
redshift, more faithfully tracing the dark matter halo mass 


function slope itself (e.g., Genel et al. 2014 Furlong et al. 


2015), while others produce a flatter slope at high redshifts 
via strong feedback (L u et ah]|2014| . As we detail below, 
a survey of the observational literature over the last several 
years finds little agreement on whether the stellar mass func- 
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tion was much steeper at z — 2 — 5 than it is at a = 0. In 
what follows, we use local observations to inform this ques¬ 
tion, relying on the large, uniform sample of star formation 
histories (SFHs) that are now available for local dwarf satel¬ 
lite galaxies ( Weisz et al.|2014a l together with the ELVIS N- 
body simulation suite of the Local Group ( Garrison-Kimmel] 
et al.|2 014) to connect local galaxy properties to the global 
galaxy population over cosmic time. 

Significant effort has been dedicated to characterizing 
the galactic stellar mass function (GSMF) and its evolu¬ 
tion. A good z ~ 0 reference point comes from |BaIdry et| 


al. (2012), who have used data from the GAMA survey to 


characterize the GSMF down to M* = 10 “Mq and find 
dfaint — —1.47. This result has been echoed by other recent 


work such as Moustakas et al. (20131 whose result agrees 
quite well with the |Baldry et al.| ( |2012| ) stellar mass func¬ 
tion, although they do not quote a value for af a i n t- 

Characterizing the mass function at higher redshift is 
also a priority. However, the results have been divergent, 
despite a significant effort. For example, |Perez-Gonzalez ~ei~| 
al. (20081 and Marchesini et al. (2009) both reported fairly 
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2 Graus et al. 


flat slopes with afaint — —1 to — 1.2 out beyond z = 2 , 


while Gonzalez et al. (20111 found steeper low-mass slopes 
af a i n t — —1.4 to —1.6 from redshift 2 = 4 — 6 , consistent 
with little evolution from the 2 ~ 0 value. Still a differ¬ 


ent conclusion was presented by Santini et al. (2012), who 


reported significant evolution, with a slope similar to the 
low-redshift value afaint — —1.45 at 2 = 1 but steepening to 
afaint — 1.85 at 2 = 3. 

Recent work has explored the stellar mass function to 
lower stellar masses, but there is still no consensus on the 


value of afaint and its evolution. Ilbert et al. (2013) re¬ 


port best-fit slopes afaint — —1.45 to —1.6 out to 2 ~ 


1.5 using UltraVISTA and Tomczak et al. (2014) use 


ZFOURGE/CANDELS and find similar best-fit values of 
afaint — —1.5 to —1.6 out to 2 ~ 3 with no obvious evo¬ 
lution. In contrast to this apparent lack of evolution, |Dun-| 


can et al. (2014) have used CANDELS and Spitzer IRAC 


data in the GOODS-S field to quantify the GSMF and find 
quite steep slopes (afaint — —1.75 to —1.9) at 2 = 4 — 5. 


Grazian et al. (20151 have done similar work in the CAN- 


DELS/UDS, GOODS-South, and HUDF fields and prefer 
slightly shallower values (afaint — —1.63) over the same red- 
shift range. Most recently, |Song et al. | (|2015|) have combined 
the CANDELS and HUDF fields with the deepest IRAC 
data to date to characterize the GSMF to an unprecedented 
depth of A/* ~ 10' M@ and report an even flatter slope 
at 2 = 4, afaint = —1.53, but see continued steepening to 
afaint — 2.45 at 2 = 8 . 

The idea that the stellar mass function might be getting 
steeper at early times is qualitatively consistent with the 
prevailing notion that cosmic reionization relies fundamen¬ 
tally on ionizing photons produced by a large number of faint 


galaxies 1 

Alvarez et al. 

[2012 

Kuhlen & Faucher-Giguere 

2012; Schultz et al.ii2014| Ro 

>ertson et al. 2015| Bouwens 


et al.|20151. Of course, a high comoving abundance of faint 


galaxies at early times requires that small dark matter halos 


are efficient at producing stars at this time (e.g., Madau et 


al. 20141, an idea that is qualitatively at odds with what 


must happen in the 2 = 0 universe. Specifically, if ACDM is 
correct, then the Milky Way should be surrounded by thou¬ 
sands of small dark matter subhalos (M v i r ~ 1O 9 M 0 ), far 


more than the number observed as dwarf galaxies (Moore 
|et al.|[l999| |Klypin et al.||1999| ). Solving the missing satel¬ 
lites problem requires the suppression of galaxy formation 
in small halos today. When (and if) small dark matter ha¬ 
los transition from being the sites of efficient star formation 
to less efficient star formation is a question that informs 
our ideas about the origin of strong feedback in small halos 
( |Jaacks et al.|[2013| |Wise et al.||2014| ). Unfortunately, the 
luminosities of interest are unobservably faint (Muv ~ —14 
at 2 > 6 ) with current telescopes, and likely will remain 
unobservable even in the era of JWST. 

An alternative way to investigate this question observa- 


tionally is to use local observations as a time machine (Weisz 
et al.|2014b Boylan-Kolchin et al.|20l4 2015 I, a possibility 


that is enabled by SFHs derived by deep color-magnitude di¬ 
agrams from HST ( Weisz et al.|2014a l. These precise SFHs 
can be used to calculate UV luminosities for each galaxy at 
specific redshifts ( Weisz et al.| 2014bl, and demonstrate that 
the progenitors of many Milky Way’s satellites (e.g., Fornax, 
Sculptor, Draco) had UV luminosities during the epoch of 
reionization that coincide well with the galaxies that seem 


to be required to maintain reionization at 2 ~ 8 (Boylan- 
Kolchin et al. |[20T~5 B15 hereafer). Moreover, by enumerat¬ 


ing these galaxies today, and tracking back their expected 
abundances using N-body simulations, [Boylan-Kolchin et alT| 
(2014 B14 hereafter) showed that the 2 = 8 UV luminosity 


function (LF) must break to a shallower slope than the ob¬ 
served a ~ —2 for galaxies fainter than Afijv ~ —14. The 
constraint was derived from the seemingly unavoidable pre¬ 
diction that there are hundreds of 2 = 0 surviving descen¬ 
dants of reionization-era atomic cooling halos in the Local 
Group. 

This paper adopts a similar approach as in B14, but 
now applied to constrain stellar mass functions at inter¬ 
mediate redshfits 2 = 2 — 5. The upper limit of this red- 
shift range is set by the age for reliable SFHs that can be 
achieved by isochrone fitting, which allows one to measure 
stellar mass older than ~ 12.5 Gyr (corresponding to 2 ~ 5) 


but with no finer precision (e.g., Weisz et al. 2014a I. The 


lower limit (2 = 2 ) is informed by our theoretical method¬ 
ology. Specifically, our method relies on studying the high- 
redshift progenitors of surviving subhalos in the simulations. 
We associate surviving subhalos with dwarf satellite galaxies 
and treat their progenitors as typical field galaxies at ear¬ 
lier times. In order to do this, we need to ensure that these 
systems were not satellites at the redshift of concern. Most 
sublialos of Milky Way size hosts in our simulations are ac¬ 
creted after 2 = 2 (see Section 2.1 below). After this time, at 
redshifts 2 < 2, the SFHs of satellite galaxies are likely af¬ 
fected by environmental quenching (|Fillingham et al.||2015 


Wetzel et al. 2015a Rocha et al. 2012), and therefore are 


poor proxies for typical galaxies in the universe prior to this 
epoch. By restricting our analysis to 2 > 2 we avoid this 
clear bias. 

The paper is organized as follows, Section 2 gives a sum¬ 
mary of the simulations and observational data used. Section 
3 presents in detail how we used the data to look for the ef¬ 
fects of galaxy suppression at high redshift and presents the 
main results of this paper. Section 4 discusses the implica¬ 
tions of this work, and provides predictions for high redshift 
stellar mass functions. The conclusions are presented in sec¬ 
tion 5, along with possible caveats and alternative explana¬ 
tions. The cosmology of the ELVIS suite and thus adopted 


for the analysis within this paper is WMAP-7 (Larson et al. 


[20TT1 < 78 = 0 . 801 , D m = 0.266, Da = 0.734, n B = 0.963, and 
h = 0.7) 


2 DATA 


2.1 Simulation Data 


We make use of the ELVIS suite (Garrison-Kimmel et al. 


2014) of N-body simulations, which is composed of 12 pairs 


of Milky Way-Andromeda analogs along with 24 mass- 
matched isolated systems for comparison. We use the 10 
best-matched Milky Way-Andromeda analogs in this work 
(see Garrison-Kimmel et al. 2014 for details). The simula¬ 
tions are complete to peak halo masses of A/haio = 6 x 10 7 
Mq, well below the mass range that can host known (classi¬ 
cal) satellite galaxies of the Milky Way (e.g. Boylan-Kolchinj 
et al.| | 2012 ), which are the points of comparison in this pa¬ 
per. Crucial to our analysis is the ability to connect bound 
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Local Group constraints on high-redshift stellar mass functions 3 


subhalos at z = 0 to their progenitors at higher redshift. 
We do this using the merger trees provided in the public 
ELVIS data basq 1 1 which were constructed using the Rock- 


star (Behroozi et al.|2013) halo finder. 

For our high redshift analysis we use only subhalos that 
survive as bound substructures at z — 0 within 300 kpc of ei¬ 
ther host. We make no assumption as to which of the pairs on 
the ELVIS catalogs would be the Milky Way or Andromeda. 
This assumption is fair because the mass estimates of the 


Milky Way and Andromeda are the same within errors (van 
der Marel et al.|2012 Boylan-Kolchin et al.|2013 |. 


One of our main assumptions in this work is that it 
is fair to treat the main progenitors of subhalos as typical 
galaxies prior to the time they first became subhalos. We 
know that classical dwarf satellite galaxies are more likely 


to be quenched than their counterparts in the field (Ma¬ 


teo 1998 and we therefore want to avoid this bias in our 
analysis. 


Using the ELVIS simulations, Wetzel et al. (2015b) 
found that ~ 84% of subhalos in the relevant mass range 
were accreted after z = 2 and ~ 98% were accreted after 
2 = 3. Our own independent analysis yields similar results 
(see also Fillingham et al. 20151. This fact motivates our use 
of 2 = 2 as the lower bound on the redshift range. We note 
that this choice is also consistent with direct estimates for 
the accretion times of the classical Milky Way satellites from 
Rocha et al. ( |2012| ), who find that only one of the Milky Way 
satellites in our sample (Ursa Minor) is consistent with hav¬ 
ing been accreted before 2 = 2, with an infall time between 
8 and 11 Gyr ago. 


2.2 Observational Data 

2.2.1 Milky Way Satellite Star Formation Histories 

The star SFHs of Milky Way satellites are taken from the 
data sets of (Weisz et al. 2014a[ W14) and Weisz et al. (2013 


W13), which rely color magnitude diagrams from archival 
HST/ WFPC2 observations. We refer the reader to |Weisz] 
et al. (2014a I and Dolphin (2002) for an in-depth discus¬ 


sion of the techniques used. In what follows we assume that 
the SFHs are applicable to the entire galaxy, even though 
in many cases the CMDs were derived from fields that do 
not cover the entire galaxy. Instead, we use the normalized 
SFHs from W14 and W13 to the 2 = 0 stellar mass of each 
Milky Way satellite using the values given in |McConnachie| 
(2012), which were calculated based on integrated light, and 


assuming a mass to light ratio of one. 

One of our major concerns is that our comparison set 
be complete. For this reason, we restrict ourselves to the 
12 dwarfs currently within 300 kpc of the Milky Way and 
with AL* =0 > 2 x 10 5 Mg: Carina, Canes Venatici I, Draco, 
Fornax, Leo I, Leo II, LMC, Sagittarius, Sculptor, Sex¬ 
tans, SMC, and Ursa Minor. The SFHs in the W13/W14 
data set include all of these bright Milky Way satellites ex¬ 
cept Sextans. In our fiducial analysis we include Sextans 
(M* =0 = 4.4 x 10 s M 0 ) by assuming its normalized SFH is 
the same as that of Carina (A/*” 0 = 3.8 x 10 5 Mg). We find 


1 http: //localgroup.ps.uci.edu/elvis/ 

2 Ultrafaint galaxies may likely be quenched, even in the field, 


that alternative choices for Sextans do not strongly alter our 
results. 


2.2.2 High redshift stellar mass functions 

We will focus our comparisons to measurements of high- 
redshift GSMF presented over the last two years. For our 
2 = 2 — 3 constraints we will normalize at the high-mass end 


to the results of Ilbert et al. (2013 113 herafter), who have 
used UltraVISTA to measure the GSMF to AL* > 10 10 Mg; 
and also to those of Tomczak et al. |2014j T14 hereafter), 
who have pushed a factor of ~ 2 deeper at the relevant 
redshifts using ZFOURGE/CANDELS. 

At 2 = 4 — 5, we will use the HST- CANDELS + 


Spitzer/ IRAC results from Duncan et al. (2014 

D14 here- 

after) and the more recent (deeper) study from 

Song et al. 

2015 S15 hereafter), who reach AL* ~ 10 7 IV 

q (vs. 10° 

VI 0 for D14) at these redshifts. |Grazian et al. (2015| have 


presented a similar analysis and their fits are intermediate 
between those of D14 and S15; we focus on the latter two in 
our comparison in order to bracket the most recent results. 


3 MOTIVATION AND METHODOLOGY 

Figure 1 provides an overview of the strategy of this work. 
The black dashed line in the upper-left panel shows the av¬ 
erage mass function of bound halos within 300 kpc of the 
ELVIS hosts at 2 = 0 (hereafter called subhalos). The red 
line shows the mass function of each bound subhalo’s main 
progenitor at 2 = 4. The shaded region shows the full scatter 
in progenitor mass functions over all 20 ELVIS hosts. We see 
that each Milky Way host should have ~ 200 bound halos 
within 300 kpc that had progenitors more massive than 10 s 
Mg at 2 = 4. Given that the number of satellite galaxies 
of the Milky Way with AL* > 10 5 Mg is significantly lower 
than this (12), we know that any GSMF that requires 10 s 
Mg halos to host A/* > 10 5 Mg galaxies at 2 = 4 would be 
inconsistent with known properties of the Milky Way. The 
aim of this work is to use this basic tension to inform our 
understanding of the on the low-mass end of the GSMF at 
high redshifts, place limits on the asymptotic slope afamt- 

The upper right panel of Figure 1 compares the best- 
fit GSMFs of D14 (with afaint = —1.89) and S15 (with 
a faint = —1.53) at 2 = 4. Table 1 lists the preferred 
Schechter parameters from D14 and S15; the lines in Fig¬ 
ure 1 show these fits and become dashed beyond the last 
published data point. 

Note that the comoving number density of AL* = 10 s 
Mg galaxies in the D14 extrapolation exceeds 10 Mpc -3 . 
We can infer right away that such an extrapolation is prob¬ 
lematic, as only halos with virial masses below 10 8 Mg at 


owing to reionization (Bullock et al. 2000 Ricotti &; Gnedin 2005 I 


these epochs have number densities this high (e.g., Lukic et 
[ah 2007). That is, the only way to have a galaxy population 
this abundant at 2 = 4 is to associate them with halos of 
mass < 10 8 Mg. As we have just discussed, the 2 = 0 de¬ 
scendants of such galaxies would drastically over-populate 
the Milky Way’s local volume. 

To add some precision to this constraint we rely on 
abundance matching, and specifically assume that there is a 
one-to-one relationship between the most massive halos and 
the most massive galaxies at 2 = 4 (we discuss potential 
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M ha io (M 0 ) M* formed before z = 4 


Figure 1 . Top left: Cumulative halo mass functions for all bound objects within 300 kpc of simulated Milky Way analogs from ELVIS. 
The black dashed line shows the average mass function at z = 0. The red line shows the average mass function for main branch halos 
at 2 = 4 that survive as independent halos to the present day. The shaded band corresponds to the full scatter over all 20 systems. Top 
right : GSMFs from D14 and S14, along with their published Schechter fits extrapolated down to the stellar masses of relevance for this 
paper. Bottom left : The stellar mass vs. halo mass relations that result from the same pair of GSMFs at z = 4. The shaded grey band 
in the lower left panel corresponds to the halo mass scale where galaxy formation may be suppressed by the UV background (|Okamoto| 
|et al.||2008| . Bottom right : The implied cumulative stellar mass functions before z = 4 that result from combining the upper left and 
bottom left panels. The data points are measurements for Milky Way satellites from W13 and W14 with Poisson error bars. We see that 
an extrapolation of the D14 GSMF would drastically over-produce the satellite population of the Milky Way. The only way the D14 
normalization can be compatible is if the GSMF truncates sharply below M* ~ 10 7 Mq. The S15 extrapolation, with its flatter faint-end 
slope, is remarkably consistent with satellite counts around the Milky Way. 
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M* (M 0 ) formed before z — 2 M* (M 0 ) formed before z = 5 

Figure 2. Projections of the stellar mass function at low masses from z = 2 to z = 5 plotted along with present day measurements of the 
stellar mass function from various sources. The solid lines are best fit faint end slopes derived by fixing the other Schechter parameters 
to literature values and performing a chi-squared minimization. In general the best fit faint end slope does not differ much from those in 
the literature, and the faint end slope appears to be consistently near -1.5 for our redshift range. 


caveats with this assumption in §5). The lower-left panel of 
Figure 1 shows the implied halo mass / stellar mass relation 
for both GSMFs. We have used t hglSheth et al. (20011 mass 
function for halo abundances herefj These relations allow us 
to assign a stellar mass to each progenitor subhalo mass at 
z = 4 and then to construct a prediction for satellite stellar 
mass functions restricted to stars formed before z — 4. 

The culmination of this exercise is shown in the bottom 
right panel of Figure 1. The black points show the cumula¬ 
tive function of galaxies within 300 kpc of the Milky Way 
as a function of M* =4 - their mass in ancient stars formed 
prior to z = 4 (from W13, W14). We note that the overall 
count is low (10 with M* =4 > 10 5 M@) and that the mass 
function is fairly flat N oc (M* =4 ) 4 ^ 3 . These counts are in 
stark contrast to what would be expected if an extrapolation 
of the D14 GSMF were to hold at z = 4 (blue), as this would 
give a steeply rising mass function (N oc M* =4 ) and result 
in > 400 galaxies with M* =4 > 10 s Mg. The shading cor¬ 
responds to the full halo-to-halo scatter in our simulations 
and the solid line shows the median. If instead we assume 
the S15 fit holds (magenta), the result is much more in line 
with the Milky Way data. Another possibility (not shown) 
is that the D14 normalization is correct but that the stellar 
mass function at z = 4 sharply breaks at M* < 10 7 Mg 
so that the number density at M* ~ 10 ’ Mg is something 
closer to the S15 extrapolation in the upper right panel. 


4 CONSTRAINTS 

Figure 2 presents the results of a similar exercise to that 
performed in the previous section but now at z = 2 (left) 


and rely on the hmf Python package I Murray et al.|2013 i. 


and z = 5 (right). As in the lower right panel of Figure 1, the 
black points show the cumulative stellar mass functions for 
Milky Way satellites within 300 kpc, specifically the stellar 
mass formed prior to z = 2 (left, M* =2 ) and z = 5 (right, 

Mr 5 ). 

The dashed lines in the left panel of Figure 2 shows the 
implied (median) cumulative mass function, N(> M* =2 ), 
for the fiducial GSMF at z = 2 from 113 (green dash) and 
T15 (red dash). Note that neither of these data sets were 
able to reach down to the low stellar masses shown in the 
plot; the results presented correspond to extrapolations of 
their double-Schechter fits (with parameters listed in Table 
1). Note that 113, motivated by their fits at lower redshift, 
have fixed Qf a int = —1.6 in their fiducial GSMFs. If this 
were to hold true, the left-hand side of Figure 2 shows that 
this would predict typically ~ 55 galaxies within 300 kpc of 
the Milky Way with more than 10 s M@ in stars older than 
about 10.4 Gyr (z = 2 formation). This is about a factor 
of five more than observed, which is larger than the halo- 
to-halo scatter expected from the ELVIS suite (shaded red 
band, see below). The T14 line, on the other hand, agrees 
fairly well with the Milky Way data, owing primarily to its 
flatter slope, af a i n t = —1.41. 

The red solid line and associated shaded band in Figure 
2 shows the preferred fit for af a j n t (—1.42) that results from 
fixing the other parameters of the T14 GSMF and then fit¬ 
ting the implied bound descendants to the cumulative M* =2 
data from Milky Way satellites. The result is almost indis¬ 
tinguishable from the mass function measured by T14 at 
z = 2. This is remarkable given that we are performing a 
consistency check ~ 10 Gyr after the epoch observed and 
~ 3 orders of magnitude lower in stellar mass than the limit 
of their survey. If instead we do the fit to the differential 
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Table 1. Schechter parameters for the reference stellar mass functions used in this work, along with the calculated best fit faint end 
slopes using local group data. Literature stellar mass functions from Ilbert+ 13 and Tomczak+ 14 were fit using a double Schechter 
function. Note that at z = 2 the literature fits for Ilbert+ 13 and Tomczak+14 are actually two separate bins (1.5 < z < 2.0 and 2.0 
< z < 2.5) which we simply averaged to get the parameters listed here. The same procedure was used at 2 = 3 for the Ilbert+ 13 data 
(using bins 2.5 < z < 3.0 and 3.0 < z < 3.5) and for Tomczak+ 14 we used the last data bin (2.5 < z < 3.0). Rows 2 and 5 through 8 
are the literature values of the Schechter fit, while Rows 3 and 4 represent the best fit slopes using local group data and keeping the rest 
of the Schechter fit paramters fixed to their literature values. 


Normalized to 113 (UltraVISTA DR 1) 


z 
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log(M*) 

4 >* 
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2 
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Table 2. Similar to Table 1 except at z = 4 and 5. In this case, both literature reference use a single Schechter fit in place of the double 
Schechter fit used in the Table 1 references. Our best fit slopes were calculated in the same manner, by keeping M* and 4'+ fixed and 
leaving a as a free parameter. 


Normalized to S15 (CANDELS/HUDF) 

z 
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log (M*) 
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4 
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1,DU -0.05 

10.68 
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mass function we find a preferred value afaint = —1.37 for 
the T14 normalization. 

The right panel of Figure 2 shifts the focus to z = 5, 
where the dashed lines show the implied Milky Way mass 
functions in stars formed prior to this redshift assuming the 
D14 (blue dashed) and S15 (magenta dashed) GSMFs. The 
D14 fit, with a steep slope afaint = —1.74, over-predicts the 
known count by an order of magnitude. The S15 GSMF at 
z = 5 is much more consistent with what we see locally. The 
solid magenta line and associated shading (full halo-to-halo 
scatter) shows the resultant best-fit afaint (= —1.62) when 
normalized to the other Schechter parameters in S15 (see 
Table 1). The result is similar to the preferred value in S15 
(—1.67), though slightly flatter. Given the significant host- 
to-host scatter, the difference between our best-fit afaint and 
the S15 value is not significant, as we now discuss. 

Table 1 lists our best-fit afaint values when normalized 
to 113 and T14 at z = 2 and z = 3 for both cumulative 
mass function fits (column 2) and differential mass function 
fits (column 3). Table 2 lists our best-fit afaint values when 
normalized to D14 and S15 at z = 4 and z = 5. The differ¬ 
ential fits weight the existence of the (likely rare) Magellanic 
Clouds less than do the cumulative fits and therefore typi¬ 
cally prefer slightly flatter slopes (at the 0.03 — 0.05 level). 
For all of the fits we perform we use the average of the cu¬ 
mulative stellar mass functions derived from ELVIS to com¬ 


pare to the W14 sample in order to compute the best fit 
faint end slope. We also compute errors on the faint end 
slope by fitting to the upper and lower 68 % range of the 
ELVIS hosts (e.g., consider the upper and lower ranges of the 
red band plotted in the upper left panel of Figure 1). This 
provides an accounting of the expected cosmic variance in 
sublialo/satellite counts from halo-to-halo. 

Figure 3 summarizes our constraints on afaint as a func¬ 
tion of redshift and compares our values to recent measure¬ 
ments reported in the literature. We specifically plot the av¬ 
erage of the four best-fit afaint values listed in Tables 1 and 2, 
with errors that treat variance among those four preferred 
values as systematic combined with the halo-to-halo vari¬ 
ance errors listed for each of the four fits. The specific values 
we prefer are: a fai nt = (-1.42±°X -1-531^5, —1.52±g'gl, 
-1.571°;°®) at z = (2, 3, 4, 5). 


5 DISCUSSION 

In this paper we have used near-field observations (deep 
CMD-based SFHs from W13 and W14) together with 
ACDM simulations to inform our knowledge of the galaxy 
stellar mass function (GSMF) at early times (z = 2 — 5). The 
power of this approach is that it utilizes a huge lever arm 
in stellar mass, exploring a mass regime two orders of mag¬ 
nitude below what is currently possible via direct counts. 
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Figure 3. Our preferred fits for af a mt as a function of redshift (black points) along with recent measurements of af a ; nt from various 
direct-imaging surveys. The redshifts of the |Song et al. (2015] andjGrazian et al. ( |2015| > results have been artificially offset for clarity. 
We include only one point from Ilbert et al. (2013) at z ~ 2 because the faint end slope they adopt is fixed to -1.6 based of their results 
at 2 = 1.5. Our black points were derived by simply averaging the best fit slopes from Tables 1 and 2 at each redshift. Our error bars 
include the variance among those four slopes as well as an accounting of halo-to-halo scatter in subhalo abundances. 


Our results suggest that the low-mass slope of the GSMF 
must remain relatively flat with only slight steepening over 
the redshift range 2 = 2 — 5, from — — 1.42+q'qj to 

~ — at 2 = 5. These constraints are strongly 

at odds with some published GSMFs, including published 
estimates that would over-produce local galaxy counts by 
an order of magnitude or more (see, e.g., the lower right 
panel of Figure 1). 

While the discriminating capacity of this approach is 
potentially very powerful, there are a few caveats that we 
must consider. The first is the concern that our simulations 
may not provide a fair representation of the Milky Way and 
its local environment. The ELVIS hosts we use have virial 
masses that vary over the range Mh a i 0 = 1 — 2.4 x 10 12 
Mq. This spans the range of most recent measurements but 
there are some estimates that fall as low as 7 x 10 11 Mq 
(see, e.g. |Boylan-Kolchin et al.|[2013| |Wang et al.||2015] . If 
the Milky Way halo mass is really this low, then we would 
expect roughly a factor of ~ 2 fewer progenitors at fixed 
halo mass within 300 kpc of the Milky Way. This would shift 
our median predictions for bound satellite counts down by 
roughly the same factor for any given GSMF, which is within 
the shaded regions shown, e.g., in Figure 2. Our preferred 
faint-end slopes would similarly be systematically steeper, 
roughly afaint —> ctfaint — 0.05. Such a shift would not change 
our conclusion that the preferred low-mass behavior of the 
113 and D14 GSMFs are too steep to be compatible with 
local counts and would not change our conclusions that af a i n t 
is evolving weakly to 2 = 5. 


A similar concern is associated with the fact that our 
simulations are pure N-body and lack the additional gravi¬ 
tational potential associated with a central disk. The extra 
tidal forces associated with the disk will likely liberate mass 
from some subhalos such that our N-body runs over-predict 


bound subhalo counts (D’Onghia et al. 2010 Brooks & Zolo¬ 


tov] 20141. We suspect that this effect will marginally affect 


our work because we count bound structures based on their 
progenitor masses prior to accretion, not on their masses at 
2 = 0 (which will certainly be lower than predicted from 
N-body for subhalos that orbit close enough to the disk). 
Progenitor masses will be largely unaffected by the presence 
of a disk. A more important concern is that more satellites 
will be completely destroyed when the disk potential is in¬ 
cluded, thus liberating their stars to the stellar halo (e.g., 
Bullock fe Johnston|2005 1 . Given that the total stellar mass 


of the Milky Way’s stellar halo is quite modest (~ 10 J Mq) 
compared to the integrated mass of the stellar mass func¬ 
tions shown in, e.g., Figure 2, there is not much room to 
hide additional mass in tidally destroyed galaxies that are 
not already accounted for in current N-body based models 
( Lowing et al.| [2015). More work needs to be done in order 
to solidify this expectation (see, e.g. Brooks et al.|[2013 for 
a promising approach). 


Another issue that could affect our results is that we 
have chosen to use one-to-one abundance matching at 2 = 
2 — 5 in order to assign ancient stellar masses to the progen¬ 
itors of our bound subhalos at 2 = 0. Imposing significant 
scatter in the M* — Mh a io relation could drastically affect our 
constraints. By definition, any imposed scatter would result 
in a GSMF that is identical (on average) to those used in 
the paper; once scatter is imposed, the median M* — Mhaio 
needs to become steeper in order to result in the same over¬ 
all number density (e.g., Garrison-Kimmel et al, in prepa¬ 
ration). This means that the median constraints on ataint 
would remain the same, but the uncertainty caused by host- 
to-host scatter would increase accordingly. In order to better 
understand the effect scatter would have on our results we 
performed a simple experiment, calculating cumulative stel¬ 
lar mass functions in the same manner, but allowing the 
stellar mass given by our abundance matching relations to 
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scatter by 0.2 dex according to a log-normal distribution. 
The resulting best fit faint end slope shifts as anticipated, 
however the shift is small (about ± 0.01) so we believe our 
results are robust unless the scatter in abundance matching 
is much larger than 0.2 dex. 

While a number potential caveats need to be under¬ 
stood before we can gain further accuracy and precision 
on our constraints, this paper (along with B14 and B15) 
has illustrated the potential power of near-field data sets in 
constraining science that is traditionally relegated to deep 
field approaches. None of the concerns we have discussed 
can account for the discrepancy between Milky Way satellite 
counts and the D14 GSMF at 2 = 4 (see Figure 1); a mass 
function as steep as afaint — —1.9 at 2 = 4 is effectively 
impossible to reconcile with our current understanding of 
structure formation. 

There are other lines of evidence from the near-field that 
can inform our understanding of the relationship between 
halo mass and stellar mass (and thus galaxy abundances) 
at high redshift. For example, it has been known for some 
time that reionization will likely shut down galaxy forma¬ 
tion below a critical halo mass scale of about ~ 10 9 Mq (Ef- 
stathiou|l992, Shapiro et al.|1994 1 and that this might leave 


imprints on the count and star-formation histories of Local 


Group galaxies that reside within those small halos (Bullock 


et al.|2000||Ricotti fe Gnedin|2 005; Gne din fe Kravtsov|2006[ 


Bovill fe Ricotti||2009 1. The shaded grey band in the lower 

left panel of Figure 1 corresponds to the halo mass scale 
where galaxy formation is expected to be suppressed by the 
external UV ionizing flux at 2 = 4 according to |Okamoto| 
et al. (2008). The D14 abundance-matching relation would 


suggest that we should begin seeing the effects of reioniza¬ 
tion suppression in galaxies as massive as M* = 10° Mg at 
2 = 4; however, no systematic truncation in star formation 
at this epoch is seen in the CMD-derived SFHs of galax¬ 
ies this massive (W13, W14). The steeper S15 relation, on 
the other hand, would suggest that UV suppression sets in at 
the stellar mass regime of ultra-faint dwarfs (M* ~ 10 4 Mg). 
This expectation, which only holds for flatter stellar mass 
functions as preferred by S15 (and our own work), agrees 
very well with recent observations that have seen uniformly 
old stellar populations in ultra-faint galaxies ([Brown et al. 


2014|) and recent simulations (|Onorbe et al.|[2015| | Wheeler 


et al. 2015 I as well. More data will be needed to confirm 


whether this dividing line between uniformly ancient stel¬ 
lar populations and continued star formation is sharp at the 
stellar mass scale of ultra-faint dwarfs, but it is clear that 
the SFHs of local dwarf galaxies have a lot to tell us about 
the physics of galaxy formation at early times. 

Finally, it is worth stressing the major assumption of 
this work, and the general concept of using the Local Group 
as a time machine, which is that the galaxies of the Lo¬ 
cal Group filled a much larger volume at high redshift, and 
thus are a representative sample of the dwarf population as 
a whole at high redshift. It is possible that this assumption 
is not correct and that dwarf galaxies that are destined to 
fall into the Local Group are somehow biased with respect 
to the bulk of the population at high redshift. Future direct 
surveys of the high-redshift universe that allow deeper con¬ 
straints on af a i n t will provide a test of this idea. If future 
counts measure different values of afaint than we anticipate 
from near-field studies, then this would help us understand 


how (and why) Local Group dwarfs would be biased relative 
to the general population of dwarf galaxies - a result that 
would be important (if true). 

The approach adopted here will be further strengthened 
by more complete observational surveys of the Local Group. 
In the near term, complete SFHs from M31 dwarfs will be in¬ 
valuable. Andromeda is believed to have a similar halo mass 
to the Milky Way, so including an additional constraint from 
this system will help us account for halo-to-halo scatter in 
subhalo populations, which dominates the error we report 
on afaint- More generally, with surveys like LSST poised to 
come on line, the next decade will almost certainly see the 
discovery of many more dwarf galaxies throughout the Lo¬ 
cal Volume, pushing our completeness limits to lower stellar 
masses and more distant radii ( Tollerud et al.|2008 l. These 
discoveries, and associated follow-up, will further lengthen 
the lever arm on stellar mass functions and increase the vol¬ 
umes within which these constraints can be applied, thus 
decreasing cosmic variance uncertainties. The same instru¬ 
ments that promise to reveal the high redshift universe di¬ 
rectly (e.g., JWST, TMT, GMT, E-ELT) will also facilitate 
indirect constraints on these epochs via archeological studies 
in the very local universe. 
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